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ABSTRACT 



Akima cubic spline and spatial Fourier series (SFS) techniques for modeling ocean 
bottom contours from bathymetric data were comparatively analyzed. SFS methods en- 
countered difficulty in ocean bottom reconstruction despite several enhancements to the 
fundamental technique. Both Akima cubic spline and SFS approaches proved unsatis- 
factory’ in reproducing first and second-order derivatives for several arbitrarily shaped 
ocean bottom contours. 

A simple reflection angle algorithm for arbitrary one-dimensional bottom models 
was examined and found to be accurate. A graphical demonstration of acoustic ray 
interaction with the ocean bottom using a variety of mathematical ocean bottom models 
and the reflection angle algorithm was conducted. 
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1. INTRODUCTION 



The primary endeavour of this thesis research was the comparative evaluation of 
existing numerical cubic spline techniques with spatial Fourier series (SFS) methods for 
the interpolative modeling of real, discrete, ocean bathymetric data. The Akima cubic 
spline was chosen for its reported ability to produce ocean bottom contour curves which 
pass through the original data while minimizing artificial oscillations in the interpolation. 
The Akima cubic spline logic was implemented via the 1MSL Version 10 double preci- 
sion computer program DCSAKM. First and second-order derivatives of the resulting 
spline were calculated from a second 1MSL Version 10 double precision program 
DCSDFR. The SFS technique was developed as a separate generic computer program 
subroutine. 

Seven differentiable mathematical bathymetric data models were employed in deter- 
mining the accuracy and relative worth of both the Akima cubic spline and SI'S repres- 
entations of the ocean bottom contour, including the first and second-order derivatives 
of these contours. A series of numerical enhancements were also applied to both mod- 
eling techniques to optimize their utility in accurately reconstructing sampled ocean 
bottoms. 

The second objective of this research was the development and evaluation of a re- 
flection angle algorithm, derived from the principles of the calculus as outlined by 
Goodman [Ref. 1], which enables the computation of reflection angles for acoustic rays 
incident upon an ocean bottom whose variation in depth is modeled as a function of 
down range. This algorithm requires exact first derivatives for angle calculations. Hence 
the apriori requirement for accurate ocean bottom modeling. 

Both of the developments outlined above were then incorporated into an existing 
acoustic ray trace program and several simple case studies investigated. 
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II. THEORETICAL ANALYSIS 



A. SPATIAL FOURIER SERIES REPRESENTATION OF OCEAN BOTTOM 
CONTOURS 

For any ocean region of horizontal range extent Z meters with a one-dimensional, 
range dependent bathymetry, the ocean bottom contour >^(z) can be modeled by the 
following spatial Fourier series (SFS): 

5 

y B (z) = | w 0 | + 2^T | w s 1 cos(2 nsf z z - Lw s ) , 0 < z < Z. (2.1) 

j=i 

The magnitude and phase angle of the complex Fourier series coefficient w s , for spatial 
harmonic 5 , are given by | wj and L\\\ , respectively. The average value for depth, which 
corresponds to zero spatial frequency or the "DC component" of the bathymetry data, 
is represented by | iv 0 | . The fundamental spatial frequency f 0 is defined as: 

/ =— 7 — ( 2 . 2 ) 

JZ o a 7 v > 

where N is the total number of bathymetric data points and T Sz is the uniform sampling 
period in the z direction in meters per sample. The highest spatial harmonic is S and is 
governed by the number of data points: 

5 = -y- , for even N (2.3) 

and 

5 = . A’ .- J - L t for odd N. (2.4) 

The first and second-order derivatives of the SFS representation of discrete ocean 
bathymetry are readily determined from Equation (2.1): 

s 

Yfl'OO = I Wj | sin(27rs/ Zo z - Lw s ) , 0 < z < Z (2.5) 

j=i 



2 



and 



2 rl 



T.S ( 2 ) ~ ~ fz. 



S 

l 



S= 1 



5 2 1 Wj | 2 cos(2 nsf z z — Lwf , 0 < z < Z. 



( 2 . 6 ) 



The complex Fourier series coefilcienls w s are determined from the spatial discrete 
Fourier transform (DFT) of the ocean bathymetric data as follows: 



N - 1 



W, 



= j~Y j y B (n T s y W™, s = -S,...,0,...,5 



n=0 



The term Uf is known as the DFT weighting factor and is given by 



where 



A'= 25+ 1 



(2.7) 



( 2 . 8 ) 



(2.9) 



when the ocean bottom sampling frequency equals the Nyquist rate. 

The reconstruction of any one-dimensional function or contour by a Fourier series 
of finite extent will be characterized by small oscilliatory deviations from the exact 
function values in the interpolated region between the discrete samples. These oscil- 
lations are often of little consequence due to the small amplitudes involved. However, 
efforts to model a simple ocean bottom profile of constant non zero slope by SFS 
methods generate significant amplitude errors at the function extremeties. This develop- 
ment highlights a particular problem that SI S methods have with apparent function 
discontinuities represented by the boundaries of the discrete data string. This anomoly, 
first documented by Willard Gibbs in 1899 and known as "Gibbs' phenomenon" (Ref. 
2], can be partially redressed by the inclusion of more data samples. This results in the 
compression of the larger amplitude oscillations toward the discontinuities. However, the 
peak amplitude of these ripples, being independent of the number of discrete data, re- 
mains constant. The problems associated with SFS methods and the "Gibbs' phenome- 
non" take on greater significance when derivatives of the bottom contour are sought. 
The ripples generated at the extremities by the Fourier series representation will be 
magnified in subsequent calculations making reliable determination of function 



3 



derivatives difficult. Strum and Kirk [Ref. 2] provide a more detailed treatment of 
"Gibbs' phenomenon". 

Lanczos [Ref. 3] argued that the impact of this phenomenon could be minimized by 
hastening the convergence of the Fourier series and proposed the following modified 
complex Fourier series coefficient tv, as a means of achieving this: 






sinfirc/S) 
sn/S 



tty ,5^0 



( 2 . 10 ) 



where vv, is the original coefficient given by Equation (2.7). 

An additional modification to Lanczos' formula was suggested by Hamming [Ref. 
4] and involved increasing the harmonic S by 1. The final expression for spatial complex 
Fourier series coefficients employing the Lanczos smoothing methodology therefore be- 
comes 



W c 



sin {snl(S + 1)) 
S7lj(S + 1 ) 



uy ,s^0. 



( 2 . 11 ) 



D. THE REFLECTION ANGLE ALGORITHM 

Acoustic rays obey the reflection laws applied in optics. Equation (2.12) relates the 
angle of incidence £ and angle of reflection V for a sound ray impinging the ocean sur- 
face or ocean bottom: 



£ = £'. ( 2 . 12 ) 

Consider a ray launched from a sound source with initial angle /f 0 and propagating 
in an isospeed ocean with a perfectly flat sea surface. When the ocean bottom is flat, the 
angle of incidence C will remain constant for successive reflections as shown in Figure 
1. The magnitude of the reflection angle is given by 

C' = /?o. (2-13) 

When the ocean bottom shoals (up slope) with constant angle of slope y, the re- 
lationship given by Equation (2.13) is invalidated and subsequent acoustic ray reflections 
become linearly dependent upon the bottom slope. This observation, reported by Officer 
[Ref. 5] and illustrated in Figure 2, is given by 

C„»(/?o-«y) ,11=1,2,3,... (2.14) 



4 



where (// = I) corresponds to the first bottom bounce with higher values for n repres- 
enting subsequent surface and bottom reflections. Negative values for indicate acous- 
tic rays which have been reflected back toward their source. Assigning a negative value 
to y in Equation (2.14) yields the relationship for £„ after an initial reflection from a 
bottom contour that deepens (down slope) with constant angle of slope y. figure 3 and 
Equation (2.15) refer to this situation: 



C* = (/?o + ">') , n= 1,2,3... (2.15) 




Figure 1. Acoustic Ray Impinging a Flat Ocean Bottom: This diagram represents 

an acoustic ray launched with initial angle of propagation in an 
isospeed ocean with a Hat ocean surface. 
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Surface 




Figure 2. Acoustic Ray Impinging a Constant Up-Slope Ocean Bottom: This dia- 

gram represents an acoustic ray launched with initial angle of propa- 
gation fi t in an isospeed ocean with a flat ocean surface. 




Figure 3. Acoustic Ray Impinging a Constant Doun-Slope Ocean Bottom: This 

diagram represents an acoustic ray launched with initial angle of propa- 
gation in an isospeed ocean with a Hat ocean surface. 
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These laws of reflection are equally valid in more general ocean models. However, 
their application is correspondingly more complex. In the case of constant sound speed 
c, acoustic rays travel in straight lines within the initial YZ launch plane making the 
calculation of ( relatively simple. However, if sound speed is made to vary with depth 
such that 



c = c(y), (2.16) 

then the acoustic ray will remain in the initial YZ launch plane but travel a curved path. 
Additionally, if the ocean bottom contour y B varies arbitrarily with down range z from 
the sound source where 



ya ( 2 - 1 7 ) 

then the relationships given by Equations (2.14) and (2.15) arc no longer relevant since 
the angle of slope y does not remain constant. 

The calculation of £ in an ocean model characterized by Equations (2.16) and (2.17) 
is achieved by representing the direction of the acoustic ray at any point along its path 
by the unit vector 

fi(y) = v(y)j + w(y) k (2. IS) 



as shown in Figure 4. 




Figure 4. Vector Representation of an Acoustic Ray: This diagram shows the 

vector representation of a sound ray travelling in an ocean medium 
where sound speed is a function of depth y. 
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The vector components v(y) and u-(y) are the direction cosines with respect to the Y and 
Z axes, respectively, describing the direction of propagation of the ray and are given by 

v(y) = cos p(y) (2.19) 

and 

w(y) = sin p(y). (2.20) 

In addition, the tangent plane and normal vector at any point P(jr 0 ,z c ) on the ocean 
bottom must be determined. 

The tangent plane at a particular point P{x^,z^) on a two-dimensional surface 
given by 



y=A*A (2.21) 

can be determined as follows (see Goodman [Ref. 1]): 

y ~ Jo =fx( x o< z o)(* ~ *o) +fz( x o» z o)( z - z o) (2-22) 



where 



fx( x o x o) 




*= Xq,Z= Zq 



(2.23) 



and 



MWq) = ~hA x ,z) I x= ^, 7= v (2.24) 

Goodman further showed the unit vector normal to the tangent plane at P(z f 0 j’ 0 ,z 0 ) to 
be 



A «± 

«l =T=“ 



(2.25) 



where 



= 



~fx( x +J ~fz( x Q’ z o) k 



(2-26) 
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and 



l /? ± I = \/Z(^o? + 1 +f^o ) 2 ■ 



(2.27) 



Substitution of the function y B given by Equation (2.17) into Equation (2.22) yields the 
following general expression for the tangent plane at any point on an ocean 

bottom which is an arbitrary function of down range z. 



^-J ; o=/b(z)(2-z 0 ). 



(2.28) 



Similarly, substitution of the function y B into Equation (2.25) provides the following 
general relationship for the vector normal to the bottom at any point P(j’ e ,z 0 ): 



where 



A 

n BL = 




A A 

"b± =+J -/b( z o) k 



and 



I n Bl I = \/ * + O-'s(’o )) 2 • 



(2.29) 



(2.30) 



(2.31) 



The angle of slope y at any point along the one-dimensional ocean floor y 5 can be 
determined from elementary vector calculus with the use of Equation (2.29) as follows: 

y = cos~ ] (y • n B1 ) (2.32) 

where y is the unit vector lying along they axis. The vector representation of an acoustic 
ray given by Equation (2. IS) can be used in defining the sound ray h t incident upon the 
ocean bottom where 



n i =v i j+w i k. (2.33) 

The angle of arrival p„ which is measured from the y axis to and is necessary in the 
calculation of the angle of incidence l, is given by 

P, = cos -1 <y .»,). (2.34) 
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Finally, with general expressions determined for y and /?„ the angle of incidence £ 
and the launch angle for the reflected ray /?, can be evaluated from the following re- 
lationships: 



£ = /?, — Ar fi y (2.35) 

and 

P r = 180° — (£ + Ay fl y) (2.36) 

where 

A y' B\ Z U -,nx 

= 7-7777-- ( 2 - 37 ) 

I y b(7 I 

Figures 5 and 6 illustrate the respective cases of an acoustic ray incident upon an ocean 
bottom with negative and positive gradients, respectively. 



10 




Figure 5. Acoustic Ray Impinging an Arbitrary Ocean Bottom: This diagram re- 

presents an acoustic ray incident upon an ocean bottom which is 
shoaling (up slope) at the point of impact. 
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deepening (down slope) at the point of impact. 
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III. OCEAN BOTTOM MODELING OF DISCRETE BATHYMETRY 



A. A COMPARATIVE ANALYSIS OF AKIM A CUBIC SPLINE AND SPATIAL 
FOURIER SERIES REPRESENTATIONS OF BATHYMETRIC DATA 

The Akima Cubic Spline (ACS) algorithm accurately models most continuous 
functions from discrete data but will frequently return inaccurate first and second-order 
derivatives. Reliable derivative calculations are crucial for reflection angle and ray 
acoustic computations and consequently an alternative modeling technique was needed. 
One of the primary focuses of this thesis, therefore, was the investigation of Spatial 
Fourier Series (SFS) methods for modeling the ocean bottom from discrete bathymetric 
data. This particular approach possessed the added attraction in that bottom roughness 
could be related to spatial frequency components. Furthermore, derivative calculations 
could be readily determined from the resultant series representation of the contour ob- 
viating the requirement to call complex external computer programs. 

Synthetic bathymetric profiles comprising 11 equidistant data points over a 1000 
meter range were generated from seven exact mathematical functions containing features 
consistent with real ocean bottoms. ACS and S I S representations of each bottom con- 
tour, including their first and second-order derivatives, were then computed. Percentage 
error calculations were made with respect to the exact function values for each model 
and its associated derivatives. When the difference between the modeled and exact data 
was of magnitude less than 1 0 6 , the percentage error was assumed to be zero. Similarly, 
percentage eriors greater than 10 6 were represented in tabular output as a series of as- 
terisks. 

The exact functions used to generate the artificial ocean bottom contours were: 

• Flat bottom (v £ (z) = 100 m ) 

• Constant downslope bottom ( y B [z) = 220.0 + 0.2z m ) 

• Constant upslope bottom ( j/ £ (z) = 420.0 — 0.2z m ) 

• One period cosine bottom ( y £ (z) = 450.0 — 250.0 cos(0.002^z) m ) 

• Modified Witch of Agnesi bottom ( y B {z) = 100.0 - — m ) 

• Half catenary bottom ( j/ £ (z) = 315.0 — 30.0(e 2 ' 600 0 + e~ 2 ' m -°) m ) 

(z - 5000. 0) 2 

• Parabolic bottom (r £ (z) = 50.0 + yr— — m ). 
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Figures 7, 8 and 9 and Tables 1 and 2 reveal shortcomings associated with the SFS 
representation of a flat ocean bottom. While both the ACS and SFS techniques repli- 
cated precisely the original contour, the latter introduced small errors in both the First 
and second-order derivatives. 

Figures 10, 11 and 12 and Tables 3 and 4 provide a comparison between the ACS 
and SFS representations of an ocean bottom which deepens (down slope) with constant 
slope. Table 3 and Figure 10 show that the ACS technique accurately modeled the ori- 
ginal function and its derivatives. However, Figure 10 and Table 4 highlight difficulties 
in the SFS representation of this same contour. This method accurately reproduced the 
original discrete bathymetric data values but introduced oscilliatorv errors in the in- 
terpolated region between them. These small oscillations around the true function values 
were propagated through the subsequent derivative calculations generating large per- 
centage errors. 

Figures 13, 14 and 15 and Tables 5 and 6 provide an analogous comparison for an 
ocean bottom which shoals (up-slope) at a constant rate. The results observed for this 
model were identical to those for the constant down slope ocean bottom. 

Figures 16, 17 and 18 and Tables 7 and 8 pertain to the ACS and SFS represent- 
ations of the one period cosine ocean bottom. Both techniques accurately reproduced 
the fundamental function shape as shown by Figure 10 but the ACS approach had 
trouble with both the first and second-order derivatives. Figures 17 and 18 and Table 7 
refer. Additionally, both the ACS and SFS reconstructions of the first and second-order 
derivatives were characterized by large errors at the contour end points. 

Figures 19, 20 and 21 and Tables 9 and 10 provide a comparison between the ACS 
and SFS representations of the modified Witch of Agnesi ocean ( see [Ref. 6] ) bottom 
model. The SFS technique was incapable of reproducing this contour or its derivatives 
as shown in the referenced Figures and Table 10. Furthermore, the oscilliatory errors 
first observed with the constant up slope and down slope ocean bottoms were again 
present. The ACS approach yielded a far more accurate reproduction of the exact bot- 
tom function. However, both the first and second-order derivative values, although bet- 
ter than those reproduced by SFS methods, were also unsatisfactory. 

Figures 22, 23 and 24 and Tables 11 and 12 relate to the ACS and SFS represent- 
ations of the half catenary ( see [Ref. 6] ) ocean bottom. The problems observed with the 
SFS reconstruction of this contour were identical to those observed for the modified 
Witch of Agnesi bottom as seen in Table 12. The ACS method again provided a true 
reproduction of the exact mathematical bottom model. In addition, the first-order 
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Figure 9. Flat Ocean Bottom : Second-Order Derivatives. 
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Figure 10. Down-Slope Ocean Bottom : Contour Reconstruction. 
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Figure 11. Down-Slope Ocean Bottom : First-Order Derivatives. 
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Figure 13. Up-Slope Ocean Bottom : Contour Reconstruction. 
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Figure 14. Up-Slope Ocean Bottom : First-Order Derivatives. 
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Figure 16 . Cosine Ocean Bottom : Contour Reconstruction. 
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Figure 19. Modified Witch Of Agnesi Ocean Bottom : Contour Reconstruction. 



27 




28 



CD 

>- 



O 

2 

o 

az 

_j 

o 

Q_ 

CC 

LJ 



CE 

CJ 

i 

CH 

LJ 



QC 

O 



CO 

LJ 

ZD 

O 



X 

CD 

LJ 





•LJj 




u|S 

-J'ac i 




olgl 


j e 

E 




g 


Sii 


S 






o 2 £ 
R B 2 



Figure 



21. Modified Witch Of Agnesi Ocean Bottom : Second-Order Derivatives. 
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Table 1. AKIMA CUBIC SPLINE RESULTS : FLAT OCEAN 
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Table 3. AKIMA CUBIC SPLINE RESULTS : CONSTANT 
DOWN-SLOPE BOTTOM 
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Table 4 . FOURIER SERIES RESULTS 
DOWN-SLOPE BOTTOM 
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Table 5. AKIMA CUBIC SPLINE RESULTS : CONSTANT 
UP-SLOPE BOTTOM. 
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Table 6. FOURIER SERIES RESULTS : CONSTANT UP-SLOPE 
BOTTOM. 
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Table 8. FOURIER SERIES RESULTS : ONE PERIOD COSINE 
OCEAN BOTTOM 
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Table 9. AKIMA CUBIC SPLINE RESULTS : WITCH OF 
AGNESI OCEAN BOTTOM 
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Table 10. rOURlER SERIES RESULTS : WELCH OF AGNESI 
OCEAN BOTTOM 
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derivative results were accurate with a uniformly small percentage error as reflected in 
Table 11. Second-order derivative values generated by the ACS approach were unsatis- 
factory, although this is not immediately evident from Figure 24 because of the plot 
scale. 

Figures 25, 26 and 27 and Tables 13 and 14 give the ACS and SFS results for the 
parabolic ocean bottom model. The ACS method matched both the original function 
and its first and second-order derivatives precisely as shown in Table 13. The SFS rep- 
resentation was reasonably accurate although some oscillations were again in evidence 
at the contour end points. These oscillations account for the large boundary errors in the 
subsequent first and second-order derivative computations as reflected in Table 14 and 
illustrated by Figures 26 and 27. 

B. NUMERICAL ENHANCEMENT OF THE SPATIAL FOURIER SERIES 
TECHNIQUE 

The results of the previous section clearly showed that the elementary SFS approach 
was unsatisfactory. Modifications to the fundamental SFS technique were therefore in- 
vestigated and applied where appropriate to improve overall performance. 

The following numerical enhancements were studied for each of the original seven 
exact mathematical models. For purposes of brevity, interim results for the modified 
Witch of Agnesi contour only have been included, as this function reflected the greatest 
improvement. 

1. Lanczos Smoothing. 

The modified expression for the SFS coefficients given by Equation (2.11) and 
incorporating Lanczos smoothing was introduced into the SFS computer algorithm and 
the modeling process repeated. 

Figure 28 and Table 15 show that Lanczos smoothing was effective in the sup- 
pression of the between data point oscillations witnessed previously in Table 10 and 
Figure 19. The corresponding reduction in the magnitude of error for the first and sec- 
ond-order derivatives can be seen in Figures 29 and 30. One undesirable effect of this 
smoothing technique was the introduction of artificial curvatures at the end points of the 
bottom contour. The cause of this phenomenon is attributed to the apparent disconti- 
nuity of the function generated by the immediate commencement and cessation of data 
at either end. 
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TECHNIQUES TOR NUMERICRL INTERPOLRTION OE Y ( R) 
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Figure 24. Half Catenary Ocean Bottom : Second-Order Derivatives. 
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Figure 25. Parabuiie Ocean Button! : Contour Reconstruction. 
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Modified Witch Of Agnesi Ocean Bottom Contour 
Reconstruction: Fourier Scries reconstruction with Lanczos smooth- 




Figure 29. 



Modified 

Derivatives: 



Witch Of Agnesi Ocean Bottom First-Order 

Fourier Series reconstruction with Lanczos smoothing. 
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Figure 30. Modified 
Derivatives: 



Witch Of Agnesi Ocean Bottom Second-Order 

Fourier Series reconstruction with Lanczos smoothing. 
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Table 11. AK1MA CUBIC SPLINE RESULTS 
CATENARY OCEAN BOTTOM 



HALF 
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Table 12. FOURIER SERIES RESULTS : HALF CATENARY 
OCEAN BOTTOM 
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Table 



. AKIMA CUBIC SPLINE RESULTS : PARABOLIC 
OCEAN BOTTOM 
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Table 14. FOURIER SERIES RESULTS : PARABOLIC OCEAN 
BOTTOM 
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Table 15. FOURIER SERIES RESULTS : LANCZOS SMOOTH- 
ING. 
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2. Lanczos Smoothing With Synthetic Data. 

The ability of Lanczos smoothing to suppress oscillations generated by the SFS 
technique gave it great utility. However, the artificial curvature at the end points of the 
contour was an undesirable effect which needed to be eliminated. It was unlikely that the 
problems associated with Lanczos smoothing could be countered by further modifying 
the expression for the SFS coefficients. Extension of the original bathymetric array by 
adding equally spaced synthetic data points was, however, viewed as a possible solution. 
The rationale of this procedure was to force the undesirable features generated by 
Lanczos smoothing beyond the actual region of interest. 

Let each data point y B . be represented as a function of incremented range such 

that 



with 



y B =y B (u) , / = o,i,2,...,(A'ZZ?/ > 7's — i) 



(3.1) 



z, = iAz , i = 0, 1 ,2,...,{NZBPTS - 1) 



(3.2) 



where A z is the distance between each data point and NZBPTS is the total number of 
original discrete data. Now let 



m s = 



i) ~ J’/?( z o) 
Az 



(3.3) 



and 



Tb( 2 (VZ/?/’7'.S-])) y B ( : (KZBPTS- 2)) .. 

m r = — . (s.4) 



Synthetic data points were then created by linear extrapolation, that is, 



y B {z_ t ) =y B {z 0 ) + m s z_ ( - , i = 1,2 ,. ..,SYNTHP (3.5) 



where 



(3.6) 



and 



yB( z (,\’ZBPTS-l+i)) — yB( z (NZBPTS-l)) + m F z {KZBPTS-\+C) < 1,2 S) i\TIIP. (3.7) 
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SYNTH P is the number of synthetic points to be generated. Figuie 31 provides a 
diagramatic representation of this methodology. 



Range of 
Interest 

0 




O = Synthetic Data Points 



0 = Original Bathymetric Data 

Figure 31. Generation Of Synthetic Data: Synthetic data points are generated 

by linear extrapolation. 

Figures 32, 33 and 3-1 and Table 16 reveal that the addition of live synthetic 
points to either end of the original data array was successful in shifting the end effects 
induced by Lanczos smoothing outside the horizontal range of interest. This resulted in 
a better curve fit for the modified Witch of Agnesi bottom and smaller errors in the 
subsequent first and second-order derivatives. One undesirable feature of the linear ex- 
trapolation method for simulating data was its tendency to suppress the natural curva- 
ture of the function being modeled. This phenomenon is observed in Figure 32. 

Generating synthetic data by linear extrapolation was a valid technique when 
the nature of the function being modeled was relatively featureless. However, it was not 
suited to functions with substantial variation in shape. Consequently, an alternative 
method involving the mirroring of the original data about the end points of the contour 
was created. This method is described by 

M ~ -/) =JT(Y) . ' = 1,2,...,SFA77/P (3.8) 
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TECHNIQUES TOR NUMERICAL INTERPOLATION OP YBIM) 




b 

§ 

2 

6 



i 



IWJ Si 



Figure 32. Modified Witch Of Agnesi Ocean Bottom Contour 
Reconstruction: Fourier Scries reconstruction with Lanczos smooth- 

ing and linear extrapolation. 
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Figure 33. Modified Witch Of Agnesi Ocean Bottom First-Order 
Derivatives: Fourier Series reconstruction with Lanczos smoothing 

and linear extrapolation. 
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Figure 34. Modified Witch Of Agnesi Ocean Bottom Second-Order 
Derivatives: Fourier Series reconstruction with Lanc/os smoothing 

and linear extrapolation. 
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Table 16. FOURIER SERIES RESULTS : LINEAR EXTRAPO- 
LATED DATA. 
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and 



•>'/i( z (,vztf/>rs-i+/)) ~ yji( z (NZDPTS-\-i)) > 1,2,. ...Si XI IIP. (3.9) 

Figure 35 illustrates this mirroring process. 



Range of 
Interest 

0 




O = Synthetic Data Points 



• = Original Bathymetnc Data 

Figure 35. Generation of Synthetic Data: Synthetic data is generated by mirror 

imaging the original data points. 



Figures 36, 37 and 38 and Table 17 give the results for the SI 'S modeling of the 
modified Witch Of Agnesi contour using synthetic data points derived from the miiror- 
ing technique. Further improvement in the overall curve lit and accuracy of the first and 
second-order derivatives was observed. However, some leveling of the bottom contour 
curvature was still prevalent. T he cause of this phenomenon was attributed to the ap- 
parent function discontinuity arising from the mirroring process. This problem is identi- 
fied in Figure 35 where mirror imaging at the right contour boundary' has created a spike 
in the data. In an elfort to compensate for this, the SFS technique has smoothed this 
feature, thereby suppressing the true function curvature at the end points of the range 
of interest. 

To resolve this problem of curvature suppression, the mirror imaging technique 
given by Equations (3.S) and (3.9) was altered as follows: 
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Figure 36. Modified Witch Of Agnesi Ocean Bottom Contour 
Reconstruction: Fourier Series reconstruction with Lanczos smooth- 

ing and mirrored data. 
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Figure 37. Modified Witch Of Agnesi Ocean Bottom First-Order 
Derivatives: Fourier Series reconstruction with Lanczos smoothing 

and mirrored data. 
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Figure 38. Modified Witch of Agnesi Ocean Bottom Second-Order 
Derivatives: Fourier Series reconstruction with Lanczos smoothing 

and mirrored data. 
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Table 17. FOURIER SERIES RESULTS : MIRROR IMAGED 
DATA. 
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yA-tb-y&q- 1)) - /- l,2,..,SYNTIIP 



(3.10) 



and 



yB( 2 (NZBPTS-\+i)) — yB( z (NZBPTS-i) * 1 ~ l’“ SYNTHP. (3-11) 

Figure 39 illustrates this modification. 



Range of 
Interest 




O = Synthetic Data Points 
0 = Original Bathymetric Data 



Figure 39. Generation of Synthetic Data: Synthetic data is generated from a 

shifted mirror imaging of the original data. 



Figures 40, 41 and 42 and Table 18 reflect only a marginal improvement in the 
modified Witch of Agnesi curve fit as a consequence of using the modified mirror imag- 
ing procedure. Additionally, for other exact mathematical models such as the one period 
cosine contour, this procedure proved detrimental, again because of curvature sup- 
pression. 

An alternative means of eliminating function discontinuities created by mirror 
imaging was investigated. This problem was most pronounced when the slope between 
successive data points increased toward the end points of the contour. A new method 
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Figure 40. Modified Witch of Agnesi Ocean Bottom Contour 
Reconstruction: Fourier Series reconstruction with Lanczos smooth- 

ing and shifted mirrored data. 
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Figure 41. Modified Witch Of Agnesi Ocean Bottom First-Order 
Derivatives: Fourier Series reconstruction with Lunczos smoothing 

and shifted mirrored data. 
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Figure 42. Modified Witch of Agnesi Ocean Bottom Second-Older 
Derivatives: Fourier Series reconstruction with Lanczos smoothing 

and shifted mirrored data. 
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Table 18 . FOURIER SERIES RESULTS : SHIFTED MIRROR 
IMAGED DATA. 
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for creating synthetic data was therefore developed to take advantage of this change in 
contour slope. Let 



yfa 1 ) 


( 3 . 12 ) 


'"10 - Az 


-yB( z \) 


( 3 . 13 ) 


~ Az 


y B ^. z ( SZBPTS — 1 )) — y u ( z ( N 7 . BPTS - 2 )) 


( 3 . 14 ) 


" lEX ~ Az 



and 



rn E2 = 



yuihKZBI'TS-2)) yii( z (NZBPrS-3)) 



Az 



(3.15) 



where again Az is the distance between each data point. If 

\m 2x I < |m 10 l 



(3.16) 



and, or 



I m E2 1 < | m El | , 



(3.17) 



then synthetic data is generated at the appropriate end of the profile as follows: 



}’b( z -i ) = >’yy(-o) + (iAz(m, 0 - /A//i» , /= 1,2,. ..,5 YNT1IP (3.18) 



and, or 



y B( z (NZBPTS-i+i)) — ysi^szBPTS- 1 )) + {iAz{m EX — iAm E )) , i — 1 ,2,...,S LA 1 IIP (3.19) 



where 



Am = am 10 — m 2x 



(3.20) 



and 



Am E = m Ex ~ m E2 . 



(3.21) 



71 



A comparison of the end points in Figures 35 and 43 demonstrate the utility of this new 
method in preserving curvature and eliminating apparent function discontinuities. 



Range of 
Interest 




O = Synthetic Data Points 



• = Original Bathymetric Data 



Figure 43. Generation of Synthetic Data: Synthetic data is generated by revers- 

ing the rate of change of slope of the original boundary data. 



Enhancing the original bathymetric data with synthetic points proved beneficial 
to the SFS modeling of the bottom contour, but no one procedure for creating this ar- 
tificial data dominated. Each of the methods investigated showed merit under specific 
circumstances. Consequently, the SFS computer algorithm was modified to incorporate 
the three methods outlined previously. If 

|A//j| = |m 10 -w 21 | <£ (3.22) 



and, or 



I A m E | = | m El - m K2 1 < £ , 



(3.23) 



where e is some arbitrarily small number, then the synthetic data was constructed by 
linear extrapolation as given in Equations (3.5) and (3.7). However, if 
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I ni 10 1 < | m 21 1 (3.24) 

and, or 

| m EX | < | m E 2 1 , (3.25) 

then mirror imaging as given by Equations (3,8) and (3.9) is used to generate the data. 
Finally, if the conditions given by Equations (3.16) or (3.17) hold then the method em- 
ploying the reversal of the rate of change of slope is activated to simulate data at the 
appropriate end point. 

C. SPATIAL FOURIER SERIES ANALYSIS OF THE AKIMA CUBIC SPLINE 
REPRESENTATION OF THE OCEAN BOTTOM CONTOUR. 

Results from the SFS technique contained significant errors in the reconstruction 
of both the original contour and first and second-order derivatives despite Lanczos 
smoothing and the optimization of methods for generating synthetic data. However, for 
the more complex ocean bottom models such as the modified Witch of Agnesi and one 
period cosine functions, SFS methods yielded the more accurate derivative values. 

The ACS method provided precise reproductions of the original test functions, but 
the subsequent first and second- order derivative values were reliable only for the flat, 
constant slope and parabolic ocean bottoms. Efforts were therefore made to incorporate 
the desired aspects of both methods. 

The ACS representation of the ocean bottom model was discretized at ten times the 
sampling frequency used to generate the original bathymetric data. Synthetic data was 
then constructed as described in Section B of Chapter 3, with the exception that A z was 
now based upon this new sampling frequency. For example, where previously 5 synthetic 
data points covered a 500 m range, now 50 synthetic points describe this same range. A 
SFS representation of the entire contour was then obtained. This approach was applied 
to each of the seven exact mathematical functions. 

Figures 44, 45 and 46 and Table 19 provide results for the SFS representation of the 
flat ocean bottom generated by the ACS method. The original contour is replicated 
precisely and the small errors observed in both the first and second-order derivatives 
approach the limits of machine precision. However, while these errors are very small, 
they are actually larger than those generated by the fundamental SFS reconstruction of 
the flat bottom without Lanczos smoothing or synthetic data. Table 2 refers. 
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Figure -4-4. Flat Ocean Bottom Contour Reconstruction: Spatial Fourier Series 

representation of the Akima Cubic Spline Hat bottom model. 
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Figure 46 . Flat Ocean Bottom Second-Order Derivatives: Spatial Fourier Series 

representation of the Akima Cubic Spline Hat bottom model. 
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Table 19. FOURIER SERIES RESULTS : ACS FLAT OCEAN 
BOTTOM 
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Figures 47, 48 and 49 and Table 20 give the results for the SFS representation of the 
constant down-slope ocean bottom generated by the ACS technique. The percentage 
errors recorded for both the bottom contour and the first and second-order derivatives 
are two orders of magnitude lower than the initial results given at Table 4. This re- 
duction, although significant, is not sufficient to eliminate residual errors likely to ad- 
versely impact on ray acoustic computations. 

Figures 50, 51 and 52 and Table 21 provide results for the SI 'S representation of the 
constant up-slope ocean bottom model. The observations made for this case are identical 
to those recorded for the constant down-slope bottom. 

Figures 53, 54 and 55 and Table 22 give results for the SFS representation of the 
one period cosine ocean bottom. The errors for both the contour reconstruction and its 
first and second-order derivatives were larger than those recorded in the earliest test 
provided in Table 8. Furthermore, the sampling of the ACS model for this bottom shape 
not only failed to assist SFS in more accurately reproducing the contour, but caused it 
to replicate anomolies in the ACS reconstruction of the first and second-order deriva- 
tives. 

Figures 56, 57 and 58 and Table 23 provide results for the SFS representation of the 
modified Witch of Agnesi ocean bottom. This contour was used to demonstrate the de- 
velopment of successive numerical enhancements to the fundamental SFS technique, and 
therefore, gives the best opportunity for a complete comparative analysis. Significant 
improvement was observed in the modeling of this function with SFS methods employ- 
ing Lanczos smoothing and synthetic data created by mirror imaging. A comparison of 
the results in Tables 17 and 23 show that further reductions in error were achieved by 
the SFS approach to modeling the ACS representation of the modified Witch of Agnesi 
bottom. However the recorded improvements in the first and second-order derivatives 
remain insufficient for the purposes of existing ray acoustics algorithms. 

Figures 59, 60 and 61 and Table 24 give results for the SFS representation of the 
half catenary ocean bottom. Significant reductions in the percentage errors of the con- 
tour were achieved. For the first and second-order derivatives reductions in percentage 
error of three orders of magnitude were common. However, the residual errors remained 
sizeable and prohibited the use of these results in further computations. 

Figures 62, 63 and 64 and Table 25 provide results for the SFS representation of the 
parabolic ocean bottom. The method of using the ACS model of the test function as 
input to the SFS procedure again proved beneficial with significant reductions in the 
magnitude of the percentage errors being observed for the contour. Reproductions of the 



78 




Figure 47. Domi-Slope Ocean Bottom Contour Reconstruction: Spatial Fourier 

Scries reconstruction of the Akima Cubic Spline down-slope bottom 
model. 
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Figure 48. Down-Slope Ocean Bottom First-Order Derivatives: Spatial Fourier 

Series representation of the Akima Cubic Spline down-slope bottom 
model. 
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Figure 49. Doun-Slope Ocean Bottom Second-Order Derivatives: Spatial Fourier 

Series representation of the Akima Cubic Spline down-slope bottom 
model. 
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TECHNIQUES TOR NUMERICAL INTERPOLATION Of YB(M) 




Figure 51). Up-Slope Ocean Bottom Contour Reconstruction: Spatial Fourier Se- 

ries representation of the Akima Cubic Spline up-slope bottom model. 
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Figure 51. Up-Slope Ocean Bottom First-Order Derivatives: Spatial Fourier Se- 

ries representation of the Akiina Cubic Spline up-slope bottom model. 
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Figure 52. Up-Slope Ocean Bottom Second-Order Derivatives: Spatial Fourier 

Series representation of the Akima Cubic Spline up-slope bottom 
model. 
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Figure 53. One Period Cosine Ocean Bottom Contour Reconstruction: Spatial 

Fourier Series representation of the Akima Cubic Spline cosine bottom 
model. 
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Figure 5-1. One Period Cosine Ocean Bottom First-Order Derivatives: Spatial 

Fourier Series representation of the Akima Cubic Spline cosine bottom 
model. 
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Figure 55. One Period Cosine Ocean Bottom Second-Order Derivatives: Spatial 

Fourier Series representation of the Akima Cubic Spline cosine bottom 
model. 
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Cubic Spline Witch of Agnesi bottom model. 
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Figure 57. Modified Witch of Agnesi Ocean Bottom First-Order 
Derivatives: Spatial Fourier Series representation of the Akiina Cubic 

Spline Witch of Agnesi bottom model. 



89 




Figure 58. Modified Witch of Agnesi Ocean Bottom Second-Order 
Derivatives: Spatial Fourier Series representation of the Akima Cubic 

Spline Witch of Agnesi bottom model. 
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Fi gure 59. Half Catenary Ocean Bottom Contour Reconstruction: Spatial Fourier 

Scries representation of the Akiina Cubic Spline half catenary bottom. 
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Figure 60. Hall' Catenary Ocean Bottom First-Order Derivatives: Spatial Fourier 

Series representation of the Akima Cubic Spline half catenary bottom. 
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Figure <31. Half Catenary Ocean Bottom Second-Order Derivatives: Spatial 

Fourier Series representation of the Akima Cubic Spline half catenary 
bottom. 
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Table 20. FOURIER SERIES RESULTS : ACS DOWN-SLOPE 
OCEAN BOTTOM 
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Table 21. FOURIER SERIES RESULTS 
OCEAN BOTTOM 



ACS UP-SLOPE 
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Table 22. FOURIER SERIES RESULTS : ACS ONE PERIOD 
COSINE OCEAN BOTTOM 
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Table 23. FOURIER SERIES RESULTS : ACS WITCH OF 
AGNESI OCEAN BOTTOM 
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Table 24. FOURIER SERIES RESULTS : 
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first and second-order derivatives were very accurate over the major portion of tiie range 
of interest, but were characterized by large errors at the outer end points of this region, 
as illustrated in Figure 64. This phenomenon is attributed to the reversal of slope method 
for generating synthetic data. 

D. ALTERNATIVE INTERIM MODELING TECHNIQUES 

Efforts to model arbitrarily shaped ocean bottom contours by SFS methods yielded 
unsatisfactory results in every case studied despite three numerical modifications de- 
signed to enhance the performance of this technique. ACS results for the contours were 
accurate, but the subsequent derivative calculations were erroneous. The need for a 
program which could accurately model both the bottom and its derivatives remained 
paramount. Without this, the precise calculation of reflection angles for acoustic rays 
impinging an undulating sea floor could not be realized. 

Motivated by the above requirements, attempts were made to salvage the ACS 
technique by utilizing its accurate contour reproduction and applying alternate numer- 
ical methods to obtain reliable values for the first and second-order derivatives. The 
following central differencing algorithms, extracted from [Ref. 7], were utilized: 

, , % ~}’b ( 2 o + 2/0 + Sr B (z 0 + h) — 8 y B (z Q — h) + y B (z n — 2h) 

y b(z o) : TTF, ( 126 ) 



and 



/'b(-o) 



-;>(~o + 2/Q + 1 (yp(z Q + h) - 3Q'»y/r n ) + 1 (H 7 ;( 2 0 - /;) -j- b (zq - 2h) 

\2h 2 



(3.27) 



where h is some arbitrarily small number, nominally 1/100 of the average spacing 
between each original data point. 

The accuracy of Equations (3.26) and (3.27) was verified in separate tests using se- 
veral mathematical functions. These relationships were then applied to the ACS recon- 
structions of the seven exact mathematical ocean bottom models listed in Section 1 of 
this chapter. The resultant first and second-order derivatives were an improvement on 
previous methods, but below the precision expected. These poor results are attributed 
to the nature of the ACS program and the manner in which the interpolating splines 
between each data point are generated. 

Having abandoned both the SFS and ACS methods for modeling the ocean bottom, 
attention was focused upon a computer routine which provided a divided differences 
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Figure 62. Parabolic Ocean Bottom Contour Reconstruction: Spatial Fourier Se- 

ries representation of the Akima Cubic Spline parabolic bottom. 
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Figure 63. Parabolic Ocean Bottom First-Order Derivatives: Spatial Fourier Se- 

ries representation of the Akima Cubic Spline parabolic bottom. 
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Figure 64. Parabolic Ocean Bottom Second-Order Derivatives: Spatial Fourier 

Series representation of the Akima Cubic Spline parabolic bottom. 
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Table 25. FOURIER SERIES RESULTS : ACS PARABOLIC 
OCEAN BOTTOM 
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interpolating polynomial fit to unordered, randomly spaced data. This program , written 
by Gerald and Wheatley [Ref. 7], was tested on the seven exact mathematical bottom 
models. Accurate curve fits were achieved; however, more care was required in sampling 
each contour to avoid artificial curvatures in the fit which are often characteristic of this 
numerical technique. The subsequent first and second-order derivatives generated by 
applying Equations (3.26) and (3.27) to the polynomial fit were precise. Therefore, these 
equations and the polynomial coefficients program formed the basis of an interim mod- 
eling technique which would enable verification of the reflection angle algorithm for an 
arbitrarily shaped ocean bottom. 
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IV. ACOUSTIC RAY REFLECTIONS FROM ARBITRARILY SHAPED 

BOTTOM CONTOURS 



The reflection angle algorithms given by Equations (2.13), (2.14), (2.15) and (2.35) 
were incorporated into an existing acoustic ray trace program. This Fortran program 
had previously been limited to the flat bottom case. Seven exact mathematical functions 
where then used to model ocean bathymetry. One acoustic ray was launched for each 
bottom contour and the computer generated angles of reflection and bottom slope 
compared with the corresponding theoretical values. An isospeed ocean was assumed in 
each case in order to interpret the results easily. 

The seven exact mathematical functions employed were: 

• Flat Bottom (y £ (z) = 1000 m) 

• Constant down-slope bottom (y B {z ) = 600.0 + 0.04z m) 

• Constant up-slope bottom (y^z) = 1000.0 — 0.04z m) 

• One period cosine bottom (y B {z ) = 250.0 — 50.0 cos(0.0002;rz) m) 

• Modified Witch of Agnesi bottom (y*(z) = 200.0 - 0 OOOU^+^OO 0 

• Half catenary bottom (j/ a (z) = 315.0 — 30.0(<? z/60UD -° + c' 2/<,ouao ) m) 

(z - 5000.0) 2 

• Parabolic bottom (y £ (z) = 250.0 4 7 J Jy ' J (j — m )- 

Figures 65 through 7S illustrate the ocean bottom models and the propagated 
acoustic ray for each of the test cases. A comparison between theory and the computer 
generated results for each model are given in Tables 26 through 32 inclusively. For each 
table, the terms y° and (° represent the exact values for the angle of bottom slope and 
acoustic ray reflection angle, respectively. The terms y° and (° represent the estimates 
generated by the ray trace program and reflection angle algorithm for these same pa- 
rameters. The bottom and surface bounces are denoted by BB and SB, respectively. 

These results were uniformly accurate with almost all of the recorded errors for the 

A 

computed reflection angles, (°, below 0.1%. Errors in the estimated bottom slope, y°, 
were also reasonably small thus verifying the central differencing methods for determin- 
ing the first-order derivatives. The parabolic ocean bottom model yielded the worst re- 
sults with errors in the reflection angle, £°, averaging 0.0S%. The modified Witch of 
Agnesi and half catenary models generated the largest errors in angle of slope, y°, but 
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the magnitudes of these angles were sufficiently small such that little impact on the 
computed angles of reflection resulted. 

Arithmetic confirmation of the reflection angle algorithm is not readily determined 
for a ray launched in an ocean medium with speed of sound varying as a function of 
depth and for an arbitrarily shaped ocean bottom. However, two such cases were inves- 
tigated as a means of visual verification of the effects of a ray incident upon an arbi- 
trarily shaped ocean bottom. The one period cosine and modified Witch of Agnesi 
models were utilized after slight modification to limit the range of depth. Figures 79 and 
80 refer. The sound-speed profile is given by Table 33. 



RECONSTRUCTED OCEAN BOTTOM BATHYMETRY 




-2.0 -1.0 0.0 1.0 2.0 3.0 4.0 



S.O 6.0 7.0 a.O 9.0 10.0 11.0 12.0 13.0 

Z t KM) 



NZBPTS: 11 



Figure 65. Flat Ocean Bottom: This profile was used to test the reflection angle 

algorithm for an isospeed ocean. 
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in an isospecd ocean with a flat bottom. 
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Figure 67. Constant Down-Slope Bottom: This profile was used to test the re- 

flection angle algorithm for an isospeed ocean. 
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Figure 68. Acoustic Ray Impinging a Constant Down-Slope Bottom: This figure 

shows an acoustic ray launched at depth 15 m with initial angle 
p 0 = 20° and propagating in an isospeed ocean with a constant down- 
slope bottom. 
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Table 26. REFLECTION ANGLES : FLAT BOTTOM WITH /?„ = 45° 



BOUNCE 


Range 

(Lin) 


ly°l 


ly° 1 


% Error 


c° 


A 

c° 


% Error 


1 BB 


0.985 


0.000 


0.000 


0.00 


45.000 


45.000 


0.00 


1 SB 


1.985 


0.000 


0.000 


0.00 


45.000 


45.000 


0.00 


2 BB 


2.985 


0.000 


0.000 


0.00 


45.000 


45.000 


0.00 


2 SB 


3.985 


0.000 


0.000 


0.00 


45.000 


45.000 


0.00 


3 BB 


4.985 


0.000 


0.000 


0.00 


45.000 


45.000 


0.00 


3 SB 


5.985 


0.000 


0.000 


0.00 


45.000 


45.000 


0.00 


4 BB 


6. 985 


0.000 


0.000 


0.00 


45.000 


45.000 


0.00 


4 SB 


7.985 


0.000 


0.000 


0.00 


45.000 


45.000 


0.00 


5 BB 


8.985 


0.000 


0.000 


0.00 


45.000 


45.000 


0.00 


5 SB 


9.985 


0.000 


0.000 


0.00 


45.000 


45.000 


0.00 



Table 27. 



REFLECT 



ION ANGLES : DOWN-SLOPE BOTTOM WITH /?„ = 20 



0 



BOUNCE 


Range 

(kin) 


1 y° 1 


ly°l 


% Error 


c° 


A 

c° 


% Error 


1 BB 


0.216 


2.291 


2.291 


0.00 


22.291 


22.291 


0.00 


1 SB 


0.494 


0.000 


0.000 


0.00 


24.581 


24.581 


0.00 


2 BB 


0.7S3 


2.291 


2.291 


0.00 


26.872 


26.872 


0.00 


2 SB 


1.136 


0.000 


0.000 


0.00 


29.102 


29.162 


0.00 


3 BB 


1.504 


2.291 


2.291 


0.00 


31.453 


31.453 


0.00 


3 SB 


1.945 


0.000 


0.000 


0.00 


33.744 


33.744 


0.00 


4 BB 


2.410 


2.291 


2.291 


0.00 


36.034 


33.034 


0.00 


4 SB 


2.961 


0.000 


0.000 


0.00 


38.325 


38.325 


0.00 


5 BB 


3.547 


2.291 


2.291 


0.00 


40.615 


40.615 


0.00 


5 SB 


4.237 


0.000 


0.000 


0.00 


42.906 


42.906 


0.00 


6 BB 


4.979 


2.291 


2.291 


0.00 


45.197 


45.197 


0.00 


6 SB 


5.851 


0.000 


0.000 


0.00 


47.487 


47.487 


0.00 


7 BB 


6.803 


2.291 


2.291 


0.00 


49.778 


49.778 


0.00 


7 SB 


7.922 


0.000 


0.000 


0.00 


52.069 


52.069 


0.00 


8 BB 


9.162 


2.291 


2.291 


0.00 


54.359 


54.395 


0.00 
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Figure 69. Constant Up-Slope Bottom: This profile was used to test the reflection 

angle algorithm 1'or an isospeed ocean. 



Ill 




/?o = 55° and propagating in an isospced ocean with a constant up-slope 
bottom. 
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Table 28. REFLECTION ANGLES : UP-SLOPE BOT1 OM WITH p 0 = 55° 



BOUNCE 


Range 

(km) 


ly°l 


ly°l 


% Error 


c ° 


A 

c 


% Error 


1 BB 


0.947 


2.291 


2.291 


0.00 


52.709 


52.709 


0.00 


1 SB 


2.476 


0.000 


0.000 


0.00 


50.419 


50.419 


0.00 


2 BB 


3.515 


2.291 


2.291 


0.00 


48.128 


48.128 


0.00 


2 SB 


4.400 


0.000 


0.000 


0.00 


45.838 


45.838 


0.00 


3 BB 


5.215 


2.291 


2.291 


0.00 


43.547 


43.547 


0.00 


3 SB 


5.909 


0.000 


0.000 


0.00 


41.256 


41.256 


0.00 


4 BB 


6.557 


2.291 


2.291 


0.00 


38.966 


38.966 


0.00 


4 SB 


7.106 


0.000 


0.000 


0.00 


36.675 


36.675 


0.00 


5 BB 


7.624 


2.291 


2.291 


0.00 


34.385 


34.385 


0.00 


5 SB 


8.059 


0.000 


0.000 


0.00 


32.094 


32.094 


0.00 


6 BB 


8.474 


2.291 


2.291 


0.00 


29.803 


29.803 


0.00 


6 SB 


S.S18 


0.000 


0.000 


0.00 


27.513 


27.513 


0.00 


7 BB 


9.149 


2.291 


2.291 


0.00 


25.222 


25.222 


0.00 


7 SB 


9.417 


0.000 


0.000 


0.00 


22.931 


22.931 


0.00 


S BB 


9,676 


2.291 


2.291 


0.00 


20.641 


20.641 


0.00 


S SB 


9.S79 


0.000 


0.000 


0.00 


IS. 350 


18.350 


0.00 
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propagating in an isospeed ocean. 
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Table 29. REFLECTION ANGLES : COSINE BOTTOM WITH /?„ = 30° 



BOUNCE 


Range 

(lun) 


lv°l 


ly°l 


% Error 


c° 


o 


% Error 


1 BB 


0.107 


0.121 


0.124 


2.48 


30.121 


30.124 


0.010 


1 SB 


0.224 


0.000 


0.000 


0.00 


30.242 


30.248 


-0.019 


2 BB 


0.341 


0.383 


0.381 


-0.52 


30.624 


30.628 


0.013 


2 SB 


0.462 


0.000 


0.000 


0.00 


31.007 


31.009 


0.006 


3 BB 


0.584 


0.646 


0.644 


-0.31 


31.652 


31.653 


0.003 


3 SB 


0.713 


0.000 


0.000 


0.00 


32.298 


32.297 


-0.003 


4 BB 


0.843 


0.910 


0.910 


0.00 


33.208 


33.207 


-0.003 


4 SB 


0.983 


0.000 


0.000 


0.00 


34.118 


34.118 


0.000 


5 BB 


1.127 


1.171 


1.173 


0.17 


35.288 


35.290 


0.006 


5 SB 


1.284 


0.000 


0.000 


0.00 


36.456 


36.463 


0.019 


6 BB 


1.446 


1.419 


1.420 


0.07 


37.878 


37.883 


0.013 


6 SB 


1.625 


0.000 


0.000 


0.00 


39.297 


39.302 


0.013 


7 BB 


1.813 


1.634 


1.632 


-1.22 


40.931 


40.935 


0.009 


7 SB 


2.023 


0.000 


0.000 


0.00 


42.566 


42.567 


0.002 


S BB 


2.245 


1.777 


1.775 


-0.11 


44.342 


44.342 


0.000 


8 SB 


2.497 


0.000 


0.000 


0.00 


46.118 


46.117 


-0.002 


9 BB 


2.766 


1.774 


1.776 


0.11 


47.893 


47.893 


0.000 


9 SB 


3.070 


0.000 


0.000 


0.00 


49.667 


49.668 


0.002 


10 BB 


3.396 


1.522 


1.522 


0.00 


51.189 


51.190 


0.002 


10 SB 


3.759 


0.000 


0.000 


0.00 


52.711 


52.713 


0.004 


11 BB 


4.144 


0.922 


0.922 


0.00 


53.633 


53.634 


0.002 


11 SB 


4.555 


0.000 


0.000 


0.00 


54.555 


54.556 


0.002 


12 BB 


4.977 


0.026 


0.026 


0.00 


54.581 


54.582 


0.002 


12 SB 


5.399 


0.000 


0.000 


0.00 


54.608 


54.608 


0.000 


13 BB 


5.812 


0.879 


0.879 


0.00 


53.728 


53.729 


0.002 


13 SB 


6.200 


0.000 


0.000 


0.00 


52.849 


52.850 


0.002 


14 BB 


6.566 


1.499 


1.498 


-0.07 


51.350 


51.352 


0.004 


14 SB 


6.896 


0.000 


0.000 


0.00 


49.851 


49.854 


0.006 
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Figure 73. Modified Witch Of Agnesi Ocean Bottom: This profile was used to 

test the reflection angle algorithm for an isospeed ocean. 
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/? 0 = 70° and propagating in an isospeed ocean with a Witch of Agnesi 
bottom. 
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Table 30. REFLECTION ANGLES : WITCH OF AGNESI BOTTOM WITH 



/?„ = 70 ’ 



BOUNCE 


Range 

(km) 


ly° 1 


1 y° 1 


% Error 


C° 


O 


% Error 


1 BB 


0.372 


0.084 


0.084 


0.00 


70.084 


70.084 


0.000 


1 SB 


0.788 


0.000 


0.000 


0.00 


70.168 


70.168 


0.000 


2 BB 


1.212 


0.248 


0.247 


-0.40 


70.416 


70.416 


0.000 


2 SB 


1.647 


0.000 


0.000 


0.00 


70.644 


70.663 


0.027 


3 BB 


2.096 


0.348 


0.348 


0.00 


71.012 


71.011 


-0.001 


3 SB 


2.563 


0.000 


0.000 


0.00 


71.359 


71.359 


0.000 


4 BB 


3.04S 


0.371 


0.372 


0.27 


71.730 


71.731 


0.001 


4 SB 


3.554 


0.000 


0.000 


0.00 


72.102 


72.102 


0.000 


5 BB 


4.0S0 


0.337 


0.336 


-0.30 


72.439 


72.438 


-0.001 


5 SB 


4.629 


0.000 


0.000 


0.00 


72.776 


72.774 


-0.003 


6 BB 


5.196 


0.275 


0.276 


0.36 


73.051 


73.049 


-0.003 


6 SB 


5.7S4 


0.000 


0.000 


0.00 


73.326 


73.325 


-0.001 


7 BB 


6.3S8 


0.211 


0.213 


0.95 


73.537 


73.538 


0.001 


7 SB 


7.009 


0.000 


0.000 


0.00 


73.749 


73.752 


0.004 


S BB 


7.644 


0.157 


0.155 


-1.27 


73.906 


73.907 


0.001 


8 SB 


8.292 


0.000 


0.000 


0.00 


74.063 


74.062 


-0.001 


9 BB 


8.950 


0.116 


0.1 15 


-0.S6 


74.179 


74.177 


-0.003 


9 SB 


9.619 


0.000 


0.000 


0.00 


74.295 


74.291 


-0.005 
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Hgure 75. Half Catenary Ocean Bottom: This profile was used to lest the re- 

flection angle algorithm for an isospeed ocean. 
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Table 31. REFLECTION ANGLES : HALF CATENARY BOTTOM WITH 



/? 0 = 60° 



BOUNCE 


Range 

(kin) 


1 y° 1 


ly°l 


% Error 


c° 


A 

c° 


% Error 


1 BB 


0.416 


0.039 


0.035 


-10.25 


59.960 


59.965 


0.008 


1 SB 


0.856 


0.000 


0.000 


0.00 


59.921 


59.931 


0.017 


2 BB 


1.294 


0.125 


0.125 


0.00 


59.800 


59.805 


0.008 


2 SB 


1.728 


0.000 


0.000 


C.00 


59.671 


59.680 


0.015 


3 BB 


2.157 


0.210 


0.214 


1.90 


59.461 


59.466 


0.008 


3 SB 


2.579 


0.000 


0.000 


0.00 


59.251 


59.253 


0.003 


4 BB 


2.995 


0.29S 


0.300 


0.67 


58.953 


58.953 


0.000 


4 SB 


3.401 


0.000 


0.000 


0.00 


58.654 


58.653 


-0.002 


5 BB 


3.799 


0.388 


0.388 


0.00 


58.267 


5S.265 


-0.003 


5 SB 


4.186 


0.000 


0.000 


0.00 


57.879 


57.877 


-0.003 


6 BB 


4.563 


0.479 


0.478 


-0.21 


57.400 


57.399 


-0.002 


6 SB 


4.926 


0.000 


0.000 


0.00 


56.922 


56.921 


-0.002 


7 BB 


5.2S0 


0.572 


0.571 


-0.17 


56.350 


56.350 


0.000 


7 SB 


5.618 


0.000 


0.000 


0.00 


56.922 


56.921 


-0.002 


S BB 


5.946 


0.665 


0.665 


0.00 


55.113 


55.114 


0.002 


S SB 


6.25S 


0.000 


0.000 


0.00 


54.447 


54.448 


0.002 


9 BB 


6.560 


0.756 


0.759 


0.39 


53.6SS 


53.689 


0.002 


9 SB 


6.845 


0.000 


0.000 


0.00 


52.929 


52.930 


0.002 


10 BB 


7.120 


0.851 


0.851 


0.00 


52.078 


52.078 


0.000 


10 SB 


7.378 


0.000 


0.000 


0.00 


51.227 


51.227 


0.000 


11 BB 


7.627 


0.941 


0.941 


0.00 


50.287 


50.2S6 


-0.002 


11 SB 


7.859 


0.000 


0.000 


0.00 


49.346 


49.345 


-0.002 


12 BB 


8.083 


1.027 


1.027 


0.00 


48.319 


41.318 


-0.002 


12 SB 


8.290 


0.000 


0.000 


0.00 


47.291 


47.291 


0.000 


13 BB 


S.490 


1.110 


1.109 


0.09 


46.182 


46.182 


0.000 


13 SB 


8.675 


0.000 


0.000 


0.00 


45.072 


45.073 


0.002 


14 BB 


8.852 


1.187 


1.187 


0.00 


43.885 


43.887 


0.005 


14 SB 


9.015 


0.000 


0.000 


0.00 


42.698 


42.700 


0.005 
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Table 32. REFLECTION ANGLES : PARABOLIC BOTTOM WITH /?„ = 70° 



BOUNCE 


Range 

(kin) 


1 y° 1 


ly° 1 


% Error 


c° 


A 

c° 


% Error 


1 BB 


1.180 


5.816 


5.798 


0.31 


64.184 


64.202 


-0.028 


I SB 


1.903 


0.000 


0.000 


0.00 


5S.368 


58.405 


0.063 


2 BB 


2.451 


3.889 


3.888 


-0.03 


54.479 


54.517 


0.069 


2 SB 


2.861 


0.000 


0.000 


0.00 


50.590 


50.628 


0.075 


3 BB 


3.217 


2.722 


2.721 


-0.04 


47.868 


47.907 


0.081 


3 SB 


3.511 


0.000 


0.000 


0.00 


45.145 


45.186 


0.091 


4 BB 


3.783 


1.859 


1.860 


0.05 


43.287 


43.326 


0.09 


4 SB 


4.021 


0.000 


0.000 


0.00 


41.428 


41.467 


0.094 


5 BB 


4.249 


1.148 


1.147 


-0.09 


40.280 


40.319 


0.097 


5 SB 


4.459 


0.000 


0.000 


0.00 


39.133 


39.172 


0.099 


6 BB 


4.664 


0.514 


0.513 


-0.19 


38.619 


38.658 


0.101 


6 SB 


4.861 


0.000 


0.000 


0.00 


38.104 


38.145 


0.108 


7 BB 


5.057 


0.088 


0.08S 


0.00 


38.192 


38.232 


0.105 


7 SB 


5.255 


0.000 


0.000 


0.00 


38.280 


38.320 


0.104 


8 BB 


5.454 


0.695 


0.694 


-0.14 


38.975 


39.014 


0.100 


8 SB 


5.665 


0.000 


0.000 


0.00 


39.670 


39.709 


0.098 


9 BB 


5.881 


1.346 


1.346 


0.00 


41.016 


41.055 


0.095 


9 SB 


6.118 


0.000 


0.000 


0.00 


42.362 


42.401 


0.092 


10 BB 


6.370 


2.092 


2.093 


0.05 


44.454 


44.494 


0.090 


10 SB 


6.661 


0.000 


0.000 


0.00 


46.546 


46.587 


0.088 


11 BB 


6.980 


3.022 


3.020 


-0.07 


49.568 


49.607 


0.078 


11 SB 


7.376 


0.000 


0.000 


0.00 


52.591 


52.627 


0.068 


12 BB 


7.844 


4.337 


4.343 


0.14 


56.928 


56.970 


0.073 


12 SB 


8.498 


0.000 


0.000 


0.00 


61.265 


61.313 


0.078 


13 BB 


9.433 


6.742 


6.669 


-1.08 


68.007 


67.982 


-0.037 
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Figure 79. Acoustic Ray Impinging A Cosine Bottom: This figure shows an 

acoustic ray launched at depth 5 m with initial angle /?„ = 85° and 
propagating in an ocean medium where sound speed varies as a func- 
tion of depth. 
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/?„ = 85° and propagating in an ocean medium where speed of sound is 
a function of depth. 
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Table 33. SOUND SPEED PROFILE 



Depth (in) 


Sound Speed (m/s) 


0.00 


1540.00 


25.00 


1524.00 


30.00 


1522.00 


36.00 


1524.00 


40.00 


1529.00 


45.00 


1532.00 


50.00 


1531.00 


60.00 


1525.75 


100.00 


1514.75 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The IMSL Version 10 Akima Cubic Spline (ACS) computer computer program had 
proven useful in modeling ocean bottom contours, but was unable to provide accurate 
first and second-order derivatives for the one period cosine and modified Witch of 
Agnesi ocean bottom models. The application of Spatial Fourier Series (SFS) methods 
as an alternative means of modeling the one-dimensional ocean floor was considered a 
viable solution to this problem. However, despite several numerical enhancements to the 
fundamental technique, including Lanczos smoothing and artificial data simulation to 
eliminate end point anomalies, Spatial Fourier Series failed to achieve the accuracy re- 
quired for subsequent acoustic ray and reflection angle computations. 

Despite the aforementioned difficulties, Spatial Fourier Series methods may yet 
prove to be a viable numerical tool in modeling ocean bathymetries, but its immediate 
utility, based upon the analysis conducted, is doubtful. 

The divided differences polynomial curve fitting Fortran routine yielded results 
which were equal to, or better than, both the ACS and SFS modeling techniques that it 
replaced. Similarly, the central differencing method for determining the first and 
second-order derivatives of the ocean bottom provided the most accurate estimates when 
used in conjunction with this alternate approach to bottom contour modeling. 
Notwithstanding the foregoing observations, this polynomial routine also had deficien- 
cies, the most notable being the introduction of artificial oscillations generated following 
an imprudent choice of discrete bathymetric data to describe the contour. 

The reflection angle algorithm developed in Chapter 2 for an arbitrarily shaped 
one-dimensional, ocean bottom were validated by the computer simulations given in 
Chapter 3. 

One of the realizations of this research was that no one numerical modeling method 
studied was definitive. Further investigation of bottom contour modeling techniques, 
such as the method of least squares, is warranted. This work represents precursory ef- 
forts in the development of a two-dimensional ocean bottom model which would enable 
more advanced study of ray acoustic theory in an ocean medium where sound-speed 
varies as a function of all three spatial coordinates. 
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